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Abstract 

We describe the 'Lie algebra of classical mechanics', modelled on the Lie algebra gen- 
erated by kinetic and potential energy of a simple mechanical system with respect 
to the canonical Poisson bracket. It is a polynomially graded Lie algebra, a class 
we introduce. We describe these Lie algebras, give an algorithm to calculate the 
dimensions c n of the homogeneous subspaces of the Lie algebra of classical mechan- 
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ics, and determine the value of its entropy limn—^ c n ■ It is 1.82542377420108 . . ., 
a fundamental constant associated to classical mechanics. 
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1 Introduction. Classes of Lie algebras. 



The class of 'simple mechanical systems' are defined by pairs (Q,V), where the config- 
uration space Q is a real Riemannian manifold and the potential energy V is a smooth 
real function on Q. The phase space T*Q has a canonical Poisson bracket and a kinetic 
energy T : T*Q — > M associated with the metric on Q. In general, the smooth functions 
on a Poisson manifold form a Lie algebra under the Poisson bracket. In the case of a 
simple mechanical system, we are given two distinguished functions, namely the kinetic 
and potential energies, and one can ask what Lie algebra they generate under the Poisson 
bracket. 

In this paper we study, not the Lie algebra generated by a particular V and T, but 
the Lie algebra defined by the whole class of simple mechanical systems. That is, one 
should think of the dimension of Q as being arbitrarily large, and the metric and potential 
energies also being arbitrary. 

This question arose out of very practical considerations of the calculations required to 
derive high-order symplectic integrators by splitting and composition, used in applications 



including molecular, celestial, and accelerator dynamics |L7| [TO]. The vector field X which 
is to be integrated is split as X = A + B, where A and B have the same properties (e.g. 
Hamiltonian) as X, but can be integrated exactly. We write exp(tX) for the time-t flow 
of X. The most common such integrator is the leap-frog method 
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<p(r) := exp(-rA) exp(rB) exp(-rA), 

where the small parameter r is the time step. 

From the Baker-Campbell-Hausdorff formula M, the map v?(t) can be represented (up 
to any power in r) as a flow exp(rX), where 

X = A + B + T 2 (lL3, [B, A}] - [A, B]])) + 0(r 4 ). (1) 

Because it is the flow of a vector field 0(r 2 )-close to the original one, the integrator 
is second order. The function X is called the modified vector field in the numerical 



integration literature [10 



For simple mechanical systems, we split the Hamiltonian as H = T + V. The flow 
of (the Hamiltonian vector field of) V can of course always be calculated easily, but 
calculating the flow of the kinetic energy T requires that Q have integrable (and even 
fairly simple) geodesies. Because the Lie algebras of Hamiltonian vector fields and of 
Hamiltonian functions are isomorphic under [Xr, Xy] = X^y,T}, there is a series formally 
identical to Eq. ([!]) involving the Hamiltonians T and V with respect to the Poisson 
bracket. 
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In the series of Eq. (|]) we see the Lie algebra generated by A and B entering. Such 
series, for example in the proof of the BCH formula, are usually considered in the context 
of the free Lie algebra L(A, B) with two generators A and B. One can in fact consider 
the more general composition 

s 

Y[ exp(a«rA) exp(6jT_B) = exp(Z) (2) 
i=i 

where Z e L(A,B). Requiring Z = t(A + B) + 0(r p+1 ) for some integer p > 1 gives 
a system of equations in the dj and 6j which must be satisfied for the method to have 
order p. In the case of general A and B, then, at each order n = 1, . . . , p there are 
dim L n (A, B) such order conditions. Here L n (A,B) is the subspace of L(A, 5) consisting 
of homogeneous elements of order n. Witt's formula states that 

dim L n (A, B) = -J2 Kd)2 n/d (3) 

d\n 

where f/,(d) is the Mobius function defined by = 1, fi(d) = (— l) fc if d is the product 
of k distinct primes, and fj,(d) = otherwise. Notice that in this case 

dim L n {A,B) ; 

n 

the dimensions grow exponentially with n. The base (2 in this case) of the exponent is 
called the entropy of L(A, B). In general, the entropy of a graded vector space 0L„ is 



limsup (dim L n ) 



l/n 



if this limit exists pl| . (We shall use generalizations of Witt's formula [[12], |19| to calculate 



the dimensions and entropies of other free Lie algebras, see Eqs. Qlq) , fllTD below.) 

In this approach it is assumed that there are no Lie identities satisfied by the vector 
fields A and B. This is reasonable if one wants the method to work for all A and B. 
However, in the case of simple mechanical systems, the Lie algebra is never free, regardless 
of T, V, or the dimension of the system. There are always identities satisfied by kinetic 
and potential energy. The simplest of these is 

{V,{V,{V,T}}} = 0. (4) 

For, working in local coordinates (q,p) with T = \p T M(q)p, and recalling the canonical 
Poisson bracket {A, B} := £\ f^g - Jg§| , we have that 

rfV 

hi 
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is of degree 1 in p, and that 

h3 

is a function of q only. So and -{V, {V", T}} commute. 

Thus, it was realized early on JIB] that in deriving high-order integrators as in Eq. (0) 
for simple mechanical systems, the order conditions corresponding to {V, {V, {V, T}}} and 
to all its higher Lie brackets can be dropped. This means that more efficient integrators 
can be designed for this class of systems. Much work has been done on this special case, 
both because of its intrinsic theoretical and practical importance, and because it allows 
such big improvements over the general case. For example, one can design special ('correc- 
tor' or 'processor') methods of the form ipip^ 1 ||, special methods for nearly-integrable 
systems such as the solar system j|, |23|, special methods involving exact evaluation of the 
forces associated with the 'modified potential' (Eq. ||) ||, and so on — see [|17] for a sur- 
vey. All of these studies rely on the structure of the Lie algebra generated by kinetic and 
potential energy. Bases for this Lie algebra have been constructed, more or less by hand, 
for small orders || ||, |20|. In particular, Murua |2(J associates a unique tree of a certain 
type to each independent order condition of symplectic Runge-Kutta-Nystrom methods 
(very closely related to the problem considered here), and enumerates these up to order 



6. (Iserles et al. [IT extend this approach to some other classes of polynomial vector 
fields.) However, a systematic description of the entire Lie algebra is clearly preferred. 

Not many classes of Lie algebras have been completely described. Here are two exam- 
ples from the literature. First, Duchamp and Krob H] completely describe all partially- 
commutative Lie algebras 

L(Ax,...,A n ; [Ai,Aj] = 0,(i,j) e C) 

where C specifies the pairs of commuting variables. Second, Kirillov, Kontsevich, and 
Molev |13| studied the Lie algebra L generated by two vector fields on 1R in general 



position, conjectured that 

y^(-l) sgn(<T) K(i),x CT (2),x <T ( 3 ),x (T (4),i/] = Vx 1 ,x 2 ,x 3 ,X4,y e L (6) 

generates all identities, and calculated the dimensions of its homogeneous subspaces and 
the asymptotic growth of their dimension. If their conjecture is true, L is a Pi-algebra 
|fj. one which the identities which hold in the Lie algebra (such as Eq. (Q)) are satisfied 
by all elements of the Lie algebra. 

Returning to the case of simple mechanical systems, it is clear that every Lie bracket of 
T and V is a homogeneous polynomial in p. Furthermore, the degrees of these polynomials 
combine in a natural way. We therefore introduce the following class ^3 of Lie algebras. 
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We use the notation [XY] := [X, Y], [XYZ] := [X, [Y, Z}], and for sets X, y, [Xy] : = 
[X,y\ :={[X,Y]:XeX,Yey}. 



Definition 1 A Lie algebra L is of class ^3 ('polynomially graded') if it is graded, i.e. 
L = @ n>0 L n , and its homogeneous subspaces L n satisfy 



Note that this implies [(L ) n+l L n ] = for all n. We call the grading of L its grading by 



For example, the Lie algebra generated by kinetic and potential energy is of class ^J3, 
where the grading is by total degree in p. The Lie algebra of all polynomial vector fields 
on a linear space is of class ^3, where the grading is by total degree. We will give more 
examples later. 

Such a grading is quite different from the natural grading of a free Lie algebra. 
Two important differences are that (i) It is not abelian. For, [L 2 , [L , L }} = while 
[L , [L ,L 2 ]] C Lq. (ii) It is not finite, in the sense that elements of L n are Lie brackets 
of unboundedly many other elements of L. For example, the bracket of any number of 
elements of degree 1 is still of degree 1. 

We also need the concept of a Lie algebra which is free in a certain class. 

Definition 2 j^j Let F be a Lie algebra of class ^3 generated by a set X. Then F is called 
a free Lie algebra in the class ^3 ; freely generated by the set X, if for any Lie algebra R 
of class ^3 ; every mapping X — > R can be extended to a unique homomorphism F — > R. 
We write F = L$p(X). 

In addition to the grading by degree, L^{X) also carries the standard grading which 
we call the grading by order, generated by order(X) = 1 for all generators I 6 I and 
order ([Y, Z]) = order (Y)+order(Z). (The term orderis chosen here because it corresponds 
to order in the sense of numerical integrators, as in Eq. ([!])). 

Because of the importance of the grading by degree for Lie algebras generated by 
kinetic and potential energy, we make the following definition. 

Definition 3 The Lie algebra L<g(A,B), free in the class ^3 ; where A has degree 2 and 
B has degree 0, is called the Lie algebra of classical mechanics. 

Two Lie algebras of class ^3 are easy to describe. First, the Lie algebra with k gener- 
ators of degree > 1 which is free in the class ^3 is just the standard free Lie algebra on k 
generators — the degrees can never decrease if the Lie algebra has no elements of degree 




'n+m— 1 



if n > or m > 0; and 



(7) 



degree. 
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0. Second, the Lie algebra with generators X = {Xi, . . . , X k } of degree and generators 
y = {Yi, . . . , Y{\ of degree 1, free in the class ^3, is y © © n >o[y n X], and only contains 
elements of degree and 1. (In both of these cases, the grading by degree is in fact 
abelian.) 

However, we want to describe the Lie algebra of classical mechanics, L<g(A, B). This 
is the simplest nontrivial it includes the essential feature of ^3 that degrees can 

both increase and decrease under Lie brackets. 

The paper is organized as follows. In Section 2, we give a construction which describes 
L<p(A, B) as the direct sum of an abelian and a free Lie algebra, both with an infinite 
number of generators. In Section 3, we enumerate the dimensions of the homogeneous (by 
order) components of L<g(A, B) and hence in Section 4 numerically compute its entropy. 
Section 5 considers special cases (e.g., of mechanical systems with Euclidean metric; these 
turn out not to be free in the class ^3) and other examples of polynomially-graded Lie 
algebras. 

2 Structure of the Lie algebra of classical mechanics. 

Let ip : L(A, B) — > £«p(A B) be the unique homomorphism from the free Lie algebra to 
the free Lie algebra of class ^3. The kernel ker if) can be thought of as the set of identities 
of L<p(A, B). For example, we showed above (Eq. ([|)) that [BBBA] G kerip. This implies 
that [CBBBA] G kerip for all C G L(A, B). However, we will see below that [BBBA] is 
not the only generator of the ideal kerip. 

Our description of L<$(A, B) is based on the following two observations. First, suppose 
one wants to describe the Lie algebra with three generators A,B,C which is free in the 
class of Lie algebras with C — 0. Since C generates all identities in this class, this Lie 
algebra is just L(A, B): one merely has to drop the generator C. To generalize this idea, 
suppose the free Lie algebra L(A, B) can be factored as @^ L(Xi) for certain generating 
sets Xi with elements in L(A, B), such that some subset y of UjXj generates all the 
identities in ^}3. Then, we have 

l<p(A B) = L(Xi\(y n Xi)) (8) 

i 

— again, we merely drop these generators. 

If y only generates some of the identities of ^3, then dropping these generators gives 
a sum of free Lie algebras which is surjectively homomorphic to L<$(A, B). This can be 
used to get upper bounds for the dimensions of the homogeneous subspaces of L<p(A, B). 

Second, given a description of L<$(A, B) as such a sum (Eq. |8|) of free Lie algebras, 
we can apply standard techniques to describe it in detail, for example to construct bases, 
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to compute its dimensions with respect to degree and/or order, and to compute the 
asymptotic growth of these dimensions. 

We begin by stating the crucial tool we shall use, the Lazard factorization of free Lie 
algebras. 

Theorem 1 \T3j] Let X and y be sets of generators. Then 

l(x u y) L{y) © L(\j n > [y n x]). 

Applying the Lazard factorization to L(A, B) with X = {A}, y = {£>}, gives 

L(A, B) = B® L(A, [BA], [BBA], [BBBA], . . .) 

where the elements [i? n A] for n > 3 are all identities in Thus, L<p(A, B) is surjectively 
homomorphic to B © L(A, [BA], [BBA]). The three generators have degrees 2 (A), 1 
([-BA]), and ([BBA]). The idea now is to eliminate this new element of degree 0. 
(Formally, the generators [-B™A], n > 3, do remain in the generating set; but they and 
all succeeding Lie brackets of them will be dropped at the final stage when we pass to 
Lrp(A,B), so we do not need to keep track of them and just indicate them by *.) This 
gives 

L(A, B) = B@ L(A, [BA], [BBA], *) 

S B © [BBA] © L(A, [BA], [BBA, A], [BBA, BA], [BBA, BBA, A]), *) 

where the generators now have degrees 2, 1, 1, 0, and respectively. Continuing in this 
way we get the following. 

Theorem 2 Let the degree of A be 2 and the degree of B be with respect to the polyno- 
mial grading (Eq. ^). Then for all k > we have the following isomorphism, 

L(A,B)^Z k ®L(A,X k ,y k ,*) 

where 

x o = 0, x k+l = x k u[y k ,A], 
y = {B}, y k+1 = [y k , x k ] u [y k , y k , a\ = [y k , x k+1 ], (9) 
z = 0, z k+ i = z k u y k , 

and * represents generators which are zero in i.e., elements of the kernel of the homo- 
morphism L(A, B) —>■ Lrp(A, B). The generating sets have the following properties: 

1. All elements ofy k and Z k have degree 0, and all elements of X k have degree 1. 

2. The Lie algebra spanned by Z k is abelian. 
3- X k = [Z k , A] . 
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4- All elements 0/34 and Z k have odd order, and all elements of X k have even order. 

5. The element of smallest order in 34 is (—l) k [[BA] k B], with order 2k + 1. 

6. The element of largest order in 34 is B k , defined recursively by B = B, B k+ i = 
[B k , [B k , A}]. It has order 2 k+l - 1. 

7. The finite sets 34 and Z k converge to infinite sets Z and X = [Z, A] in the sense 
that the sets 

{X : X G £ fc ,order(X) < n} 
are all equal for k > n/2. We have 

L(A, B) = Z @ L(A, X, *) 

and 

L^(A,B)^ Z®L(A,X). (10) 

8. The sizes of the sets X k and 34 obey the iteration 

(11) 

with initial conditions \Xq\ = 0, |3^o| = 1- This iteration generates the sequence of 
\X k \ values 

0, 1, 2, 4, 12, 108, 10476, 108625644, . . . ; (12) 
there is a constant 7 w 1.1555 such that for sufficiently large k, \X k \ = [7 2 ^]- 

Proof The iteration results from successive elimination of elements of degree 0, each 
iteration introducing only a finite number of new elements nonzero in which have 
degrees and 1. The other points then follow easily. The final description of L^(A,B), 
Eq. flTD|), follows because the generators of L(A,X, *) have degree 2 (A), 1 (X), or are 
identically zero (*). Therefore L(A, X, *) contains no elements of degree 0, so L<p(A, X) = 
L(A, X). The sequence of Eq. (|l"2"D is Sloane's sequence A001696 | 22| , which comes from 



34_|_i 


— 34 


+ 1341 


134+1 1 


= 1341 


34 + i 



the same iteration (Eq. [TTj) ; the reference there to [0] shows how to establish its doubly- 
exponential growth. • 

The rapid growth of the sets X k and 34 means that it is impossible to carry out the 
iteration exactly very far. In practice the generating set Z can be found up to any order n 
by dropping any terms of order > n as soon they appear in 34 (i-e., by quotienting all Lie 
algebras by the ideal consisting of all elements of order > n). We then have 3 ; [(n+i)/2] = 
and the iteration terminates. 
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Table 1: Elements of degree and weight < 11 (i.e., functions of q only or 'modified 
potentials' of simple mechanical systems) appearing at iteration k of Eq. (|^). The new 
elements are numbered consecutively Z±, Z 2 ,.... The degree 1 elements X n := [Z n ,A] 
also appear. 



k 


34 






order 


1 


Zi = 


B 




1 


2 


z 2 = 


[Z U X X ] 


(= \BBAX) 

V L J / 


3 


3 


Zj = 


[Z 2 ,Xi] 


(= [BBA,^]) 


5 




Za = 


\ ZjO . J\. 9 
L ^ ? ^ J 


(= [BBA, [BBA, A]]) 


7 


4 


Zc, = 


[Z 3 ,Xi] 

L ° " J- J 


(= [Z?IM, £A], BA]) 


7 




Z 6 = 


[Z3, x 2 \ 


(= [BBA, A]]) 


9 




z 7 = 


[Z3, x 3 ] 




11 




Z s = 


[Za,X 1 ] 




9 




z 9 = 


[Z4, X 2 \ 




11 


5 


Z\o = 


--[Zs,Xi] 




9 




Z n = 


--[z 5 ,x 2 ] 




11 




Z\2 = 


-- [z 6 ,x 1 ] 




11 




Zrs = 


~-[Z 8 ,X 1 ] 




11 


6 


Z\A = 


z [Zio, X\\ 
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The results of the six iterations required when n — 12 are shown in Table 1. We name 
the elements of Z Z\, Z 2 , . . . as they are successively generated by the algorithm. This 
gives a short description of the elements of (L^) n (A,B) of order < 12 in terms of 14 
elements of degree 0, 14 elements of degree 1, and 1 element of degree 2, which generate 
a total of 283 elements of weight < 12 (see Tables 2 and 3). 

3 Dimensions of the homogeneous components. 

We now turn to the enumeration of Xk and 34 by order. We introduce the generating 
functions 

00 

Xk (t) = ^|{IeI fc : order(X) = n}\t n 

n=l 

oo 

Vkit) =Y^\{Yey k : order(F) = n}\t n 

n=l 

oo 

h(t) = J2\i ZeZ k- order(Z) = n}\t n 

n=l 
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which from Eq. (ST) obey 



x 


= 0, 


m 


= 0, 




= t, 


Xk+l 


= x k + ty k 


Vk+l 


= ykXk+i, 


Zk+l 


= z k + y k . 



We can eliminate the t-dependence of this map by introducing y k = ty k and z k = tz k . 
Then z k = x k for all k and the rest of the system is 



x = 0, 

(13) 



y = t 2 



x k+ i = x k + y k , 

The polynomials x k (t) converge to a formal power series x(t). The polynomials y k (t) 
converge, again in the sense of formal power series, to 0. The power series x(t) completely 
determines the dimensions of the homogeneous components of (L<p) n (A, B) (including its 
abelian part Z, because z k (t) = x k (t)/t). We find 

x(t) =t 2 + t 4 + t e + 2t 8 + 3t 10 + 6t 12 + 12t 14 + 24t 16 + 50t 18 + 107t 20 + 232t 22 + 

508t 24 + 1124t 26 + 2513t 28 + 5665t 30 + 12858t 32 + 29 3 56t 34 + 67371t 36 + . . . ' 

(14) 

(For example, the 1 + 1 + 1 + 2 + 3 + 6 = 14 generators of weight < 12 are given in Table 
1.) Amazingly, this power series has appeared before (apparently as a curiosity) from the 
same iteration (Eq. [13]), and it appears as Sloane's sequence A045761 |^2 . 



The classical formula of Witt, Eq. (|3|), can be extended to free Lie algebras with 
more general generating sets fl2], [19|]. For any set A with generating function a(t) = 
Sn>o e "4 : order(A) = n}\t n , the dimensions c n = dimL n (^4) of the homogeneous 
components of the graded Lie algebra L(A) = {& n>0 L n (A) are given by 

Cn = ^2 ^( d ) b n/d, (15) 

d\n 

where 

-log(l-a(t)) = ^6X. 

n>0 

In Maple, one can compute the dimensions easily by c=EULERi (INVERT (a) ) (these func- 
tions are available in P2"|), where a and c are the sequences of coefficients of a(t) and c(t), 
respectively. 



10 



Table 2: Dimensions of Lie algebras graded by order. Column 2: Of the free Lie algebra 
with two generators. Column 3: Of the Lie algebra of classical mechanics, L<$(A, B) where 
A ('kinetic energy') has degree 2 in p and B ('potential energy') has degree in p, i.e. is 
a function of q only. Column 4: Number of modified potentials of order n in L<g(A,B). 
Column 5: Upper bound for maximum number of linearly independent Poisson brackets 
of order n when M = R n with the Euclidean metric, i.e. A = p T p. Column 6: As Column 
5, but V(q) is a cubic polynomial. 



n 


dim L n (A,B) 


dim(Lm) n {A, B) 


[t n+L ]x(t) 


Euclidean 


Cubic 


1 


2 


2 


1 


2 


2 


2 


1 


1 




1 


1 


3 


2 


2 


1 


2 


2 


4 


3 


2 




2 


2 


5 


6 


4 


1 


4 


3 


6 


9 


5 




5 


3 


7 


18 


10 


2 


10 


6 


8 


30 


14 




14 


6 


9 


56 


25 


3 


25 


10 


10 


99 


39 




39 


12 


11 


186 


69 


6 


69 


19 


12 


335 


110 




110 


22 


13 


630 


194 


12 


193 




14 


1161 


321 




320 




15 


2182 


557 


24 


555 




16 


4080 


941 




938 




17 


7710 


1638 


50 


1631 




18 


14532 


2798 




2787 




19 


27594 


4878 


107 


4857 




20 


52377 


8412 




8376 




21 


99858 


14692 


232 


14624 




22 


190557 


25519 




25399 




23 


364722 


44683 


508 


44460 




24 


698870 


77993 




77594 




25 


1342176 


136928 


1124 


136191 




26 


2580795 


240013 




238684 




27 


4971008 


422360 


2513 


419916 




28 


9586395 


742801 




738375 




29 


18512790 


1310121 


5665 


130199 




30 


35790267 


2310451 




2295702 




31 


69273666 


4083436 


12858 


4056416 




32 


134215680 


7218252 




7169109 




33 


260300986 


12781038 


29356 


12691109 




34 


505286415 


22638741 




22474996 




35 


981706806 


40152860 


67371 


39853452 




36 


1908866960 


71247291 




70701714 




37 


3714566310 


126559227 


155345 


125562178 




38 


7233615333 


224917313 




223099566 




39 


14096302710 


400080000 


359733 


396759314 




40 


27487764474 


711997958 




705941791 
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We apply Eq. (]T5|) to L V (A,B) = Z © L(A,X). The generating function for the 
grading by order of {A} U X is t + x(t). This gives the dimensions listed in Table 2 for 
1 < n < 40. A dramatic reduction in the dimensions compared to those of the free Lie 
algebra of rank 2 is evident. 



More generally still, Kang and Kim Jl2j consider the grading of a free Lie algebra by 
an abelian semigroup S which satisfies the finiteness condition that any s G S is a sum 
of other elements of S in only finitely many ways. Then we have 

dim L S (A) = J2^(d)b s/d (16) 

d\s 

where 

-log(l-a(t)) = ^M s 

and d\s means that there exists r G S such that dr = s, in which case we write s/d = r. 

We can use this to calculate the dimensions of L<g(A, B) with respect to the bigrading 
by order and degree. We first simplify the grading by degree, Eq. (|7p, by introducing 
degree'(x) := degree(x) — 1. Then (as long as no elements of degree enter, which now 
holds), the semigroup of the grading by degree' is isomorphic to the nonnegative integers 
under addition. Including the grading by order gives S = Z >0 x Z-°. Note that the 
finiteness condition holds for S since it holds for Z >0 . Since order(A) = 2, degree'(A) = 1, 
and degree'(X) = for all X G X, the generating function of {A} U X is ut + x(t) and we 
apply Eq. ([TED with 

bt, u = -[t n u m ] log(l -ut-x(t)). 
This gives the dimensions for n, m < 16 as shown in Table 3. 

4 Asymptotics of the dimensions and calculation of 
the entropy. 

From Eq. fli~5D , the asymptotic growth of the dimensions c n is determined by the analytic 
structure — the location and type of the singularities — of — log(l — a(t)). These correspond 
to zeros and singularities of 1 — a(t). In particular, if 1 — a(t) has a simple zero at t — a 
and no other zero with \t\ < a, then 

Cn--(-Y (17) 



n \a 



and the Lie algebra has entropy 1/a. 

The generating function of {A} UX is t + x(t). We therefore need to study the analytic 
structure of the function l — (t + x(t)). We therefore study the map of Eq. (|13|) considered 
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Table 3: Dimensions of Ltp(A, B), graded by degree m and by order n. 



n 


m 
total 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 15 


1 


2 


1 





1 


























2 


1 





1 




























3 


2 


1 





1 
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2 





1 





1 
























5 


4 


1 





2 





1 






















6 


5 





2 





2 





1 




















7 


10 


2 





4 





3 





1 


















8 


14 





4 





6 





3 





1 
















9 


25 


3 





9 





8 





4 





1 














10 


39 





9 





14 





11 





4 





1 












11 


69 


6 





20 





23 





14 





5 





1 










12 


110 





18 





37 





32 





17 





5 





1 








13 


194 


12 





46 





62 





46 





21 





6 





1 






14 


321 





42 





90 





97 





60 





25 





6 





1 




15 


557 


24 





107 





165 





144 





80 





29 





7 





1 


16 


941 





90 





229 





274 





206 





100 





34 





7 


1 


as a map 






<P 


: C 2 


^C 2 , 


(x,y) i- 


► (x + y,y( 


x + y)) 











with initial conditions x = 0, y = t 2 . If the iterates of the map converge to (x*,0) say, 
then x(t) = x*. Curiously, the map preserves the area ^dx A dy, although this plays no 
role in the analysis. 

The map (p has a line of degenerate fixed points (x,0) with eigenvalues x and 1. 
The fixed points with \x\ > 1 are unstable and one can show that the fixed points with 
\x\ < 1 are stable. The map 'remembers' its initial condition, and the function x(t) is the 
re-coordinate of the fixed point reached from initial condition (0,t 2 ). 

We can see immediately that (i) for t real and positive, x(t) is strictly increasing; and 
(ii) if the map converges then \x(t)\ < 1. For t real and positive, the sequence {y k } is 
increasing, and if there is a k such that y k > 1, then x k — > oo. Therefore we define 

(3 = M{t e R + : x k (t) -> oo}. 

Because 

Ffc+l| — \ x k\ + \Zk\, l^fc+l' = |£fc||Xfc+l|, 

the map converges in the disk {t : \t\ < (3}. 

We can get a crude bound on j3 immediately, but more detailed knowledge requires a 
numerical study of the map ip. Let t be real, let I(x, y) = y + x — 1 + \^2y, and suppose 
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x > 0, y > 0, and I(x, y) < 0. Then 

Io(p-I = y(x + y)+x + y-l + ^/2y(x + y) -(y + x-l + y/2y) 



= y(x + y) + v^fV^ - 1) 




= -V2y 3 / 2 
< 

Therefore, the orbit must stay in the bounded region x > 0, y > 0, I(x,y) < 0, with Xk 
increasing and yu decreasing. Therefore the orbit converges to some fixed point (x, 0). 
(Here the curve x = l — y/2y—y was chosen because it is a good approximation of the stable 
manifold of (1, 0).) Since 1(0, t 2 ) < for < t 2 < 2 - \/3, we have (3 > \/2-^ > 0.51. 
Better approximations of (3 can be obtained as the roots of I o <£> fc (0,/3 2 ) = (i.e., by 
requiring the fcth iterate to land in the trapping region), but these must be calculated 
numerically. On the other hand, x 2 = t 2 + t 4 > 1 if t > 0.79, so we have the bounds 

0. 51 < < 0.79. 

We have that dx(t)/dt > on [0,/3), with x(0) = and x(f3) = 1; and 1 — t is 
decreasing. Therefore 1 — t — x(t) has exactly one zero in [0,(3), and it is simple. The 
zero is a, the reciprocal of the required entropy of L<g(A,B). The numerical value of a 
can be determined by solving 1 — t — x(t) = numerically^ This gives the value of the 
entropy of Ly(A, B) as 

l/a = 1.82542377420108.... (18) 

Are there any other solutions to 1 — t — x(t) = 0? Because the coefficients of x(t) are 
all nonnegative, there can be none in the disk \t\ < a. To say more we have to proceed 
numerically. Firstly, if \xk\ and \yk\ get too large then the orbit blows up. Let 

D = {(x,y) E C 2 : \y\ > 2\x\ > 2}. 

Suppose (xk,yt) £ D. Then 

\%k+l\ > IM - \Xk\\ > \x k \ > 1 

and 

\Vk+i\ — FJfe+illl/jfel > 2|a;fc + i|, 

1. e., we have {xk+i,yk+i) ^ D. The orbit then stays in D and cannot converge — in fact, it 
must blow up doubly exponentially. The first iterate (xi,yi) = (t 2 ,t 4 ) is in D if \t\ > v2, 

1 In MATLAB, by function x = f(t); x=0; y=t~2; while y>le-16, x=x+y; y=y*x; end; x = 
1-t-x; and alpha = f solve ( ' f ' , . 5) . 
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Figure 1: Contour plot of |1 — t — x(t)\, showing its main zero at t = a = 1/1.8254. . . 
and other zeros (two sequences approaching t = (3. The unit circle is also shown. 

and the second iterate (x 2 , y 2 ) = (t 2 + 1 4 , t 4 (t 2 + 1 4 )) is in D if \t\ > 1.27202. In practice, 
if an iterate enters this region one can immediately stop the calculation and report that 
the map diverges. 

Using this criterion we computed the function x(t) numerically. See Figures 1 and 2. 
We have made the following numerical observations: 

1. The singularity of x{t) closest to the origin is at 



2. There are no zeros of 1 — t — x(t) in the disk \t\ < f3. 

3. The map converges only in a connected, simply-connected region with a fractal 



4. The function x(t) is analytic everywhere inside this region but has a square root 
singularity everywhere on its boundary. 



t = (3 = 1/1.58207912734.... 



(19) 



boundary. 



5. For each point z on the boundary, x(t) 
depending on z, as t — » z. 



1 - a(t 



— z) 1 / 2 for some constant a 



6. There is only one zero of 1 — t — x(t) in \t\ < (3. 
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Figure 2: Contour plots of \x(t)\ (left, contour interval 0.05) and argrc(t) (right). 



7. The other zeros of 1 — t — x(t) form two infinite sequences a n , a n , with Re(a n ) > f3 
for all n and lim^oo a n = f3. 

Because of the fractal nature of the boundary, we are unlikely to be able to 'solve' the 
map ip or find a in closed form. Observation (4) would imply that this boundary forms a 
natural boundary for the function x(t). Observation (5) would imply that the number of 
modified potentials of order n, [t n+1 ]x(t), is 0(n~ 3 / 2 /3~ n ). Observation (6) would imply 
that the next term in the asymptotic growth of c n = dim(L«p) n (A, B) comes from the 
square root singularity at t — (5. Indeed, by computing c n numerically for n < 80 we find 
that 

_! _ n _ J IMn-Wp-" n odd 
Cn ~ U 01 \ lMn- 3 / 2 p- n neven 

and 

[t n+1 ]x(t) ~ 0.9628n- 3/2 (3- n 
for n even. These are all consistent with the observed singularity structure of 1 — t — x(t). 
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5 Discussion 



5.1 Physical interpretation of the generators 

There is a particularly nice interpretation of L<$(A, B) = Z © L(A, X) in the specific 
case of simple mechanical systems. In local coordinates, let A = T(p) = \p T M(q)p be 
the kinetic energy, where M(q) is the inverse of the metric (or mass matrix), and and 
B = V(q) be the potential energy. The set Z consists of functions of q only, and we think 
of them as 'modified potentials'. Elements of the span of Z, 

a Z T desme{z) Z = airV + a a r 3 M (V, V) + . . . , 

and their flows, can be evaluated explicitly and used to construct high-order integrators 
of the full system T + V (see Eq. (|3|) for more terms). Now consider the generator 
X = [Z, A] G X. It is the cotangent lift of the gradient flow of the modified potential Z; 
we have X = M(q)(Z(q),p) and Hamilton's equations are 

q = M{q)Z'{q) = dw M - Hq) Z =: f(q) 

v = -f'(q) T p- 

So in a sense the modified potentials and the kinetic energy together contain a complete 
description of the Lie algebra. 

5.2 Euclidean mechanical systems 

Recall that on each manifold M, each simple mechanical system (say with kinetic energy 
T and potential energy V) generates a Lie algebra of class ^3. Therefore there is a 
homomorphism ip(M,T,V) from Lsp(A,B) onto this Lie algebra. One can ask whether 
the system (M, T, V) is in general position, i.e. if the two Lie algebras are actually 
isomorphic and ker ip(M, T, V) = 0. This is unlikely, because of the existence of identities 
such as Eq. (|6]) in Lie algebras of vector fields. One can therefore consider larger classes 
of systems and ask whether they are in general position. That is, does the class satisfy 
any identities other than those corresponding to the grading by degree, Eq. (0)? We 
conjecture that for the class of all simple mechanical systems, it does not. 

Conjecture 1 The only identities satisfied by all simple mechanical systems are those 
due to the grading by degree. That is, 

p| kerip(M,T,V) = 0. 

M,T,V 

This is best discussed by introducing a smaller class which we shall see is not in general 
position. Namely, let M = M. n with the Euclidean metric. Then in coordinates the kinetic 
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energy is T(p) = | Y^=iPi- The first few modified potentials are then 
Z 1 = V 

Z 2 = [BBA] = V'{V) 

Z 3 = [BBA, BA] = 2V"(V, V) 

Z 4 = [BBA, [BBA, A}} = 4V"(V"(V), V) 

Z b = [[BBA, BA], BA] = 2V'"{V, V, V) + W"(V"(V), V) 

where we regard the /cth derivative of V as a real-valued symmetric linear function on 
k vectors. Each modified potential of order 2n — 1 is a linear combination of the scalar 
elementary differentials of order n of V. Each such differential can be associated to a free 



tree with n nodes. (See, for example, |10 for a discussion of elementary differentials and 



trees.) The number of such trees for n > 1 is (Sloane's A000055, fl22| ) 1, 1, 1, 2, 3, 6, 11 



23, 47, 106, 235, .... This should be compared with the number of modified potentials 
in Eq. (HD, namely 1, 1, 1, 2, 3, 6, 12, 24, 50, 107, 232, .... There are three interesting 
consequences: 

(i) For n < 6, the sequences are the same. In fact, one can check that in the modified 
potentials of orders 2n — 1 < 11, all trees appear, in invertible linear combinations, 
so these modified potentials are in general independent. 

(ii) For n = 7, 8, 9, there are more modified potentials than free trees. In particular, 
only 11 of the 12 modified potentials of order 13 can be linearly independent. This 
proves that the class of Euclidean mechanical systems is not in general position. 

(iii) For n > 10, there are fewer modified potentials than free trees. In fact, the former 
have entropy 1/(3 = 1.582 . . . while the latter (since the free trees have entropy given 
by Otter's constant, 2.955. . . ) have entropy V2.955 . . . = 1.719 . . .. Thus, for large 
n, only certain combinations of the trees appear in Z. 

So far we have only considered the modified potentials Z themselves. If these are 
independent, then X = [Z, A] is independent too. However, there is still a possibility 
for extra identities to hold in the Lie algebra generated by A and X. A term of order n 
and degree m is a sum of elementary differentials of V and p, corresponding to trees with 
(n + m + 1)/2 nodes, of which m leaves are labelled p and the remaining nodes are labelled 
V. In this case we find that for (n + m + l)/2 < 7 there are always sufficient labelled 
free trees to prevent forced dependencies among the Lie brackets. For example, of the 11 
free trees with 7 nodes, there are 12, 20, 24, 18, 9, 3, and 1 trees in which m — 0, 1, 2, 
3, 4, 5, and 6 leaves are coloured p, respectively. The dimensions of the corresponding 
homogeneous subspaces of L<$(A, B) with (n + m+ l)/2 = 7 are (from Table 3) 12, 18, 20, 
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14, 8, 3, and 1, respectively. Thus, only in the case m = 0, corresponding to the modified 
potentials themselves, is a dependency forced in this way 

The algorithm for L<$(A, B), Eqs. (^) and (|T3|), can be modified to take into ac- 
count the dependencies amongst the Lie brackets in the Euclidean case. To get an upper 
bound on the dimensions and entropy of the Lie algebra in this case, we assume that 
the dependency appears only when forced. Let c n be the number of free trees with n 
nodes. At iteration k, we already have z^n := [^ 2n ]^fc elements of order 2n — 1 in Z k , and 
Vk,2n '■= \P 2n ]ijk elements of order 2n — 1 have just been created in 34- If -2fc,2n + Vk,2n > c n , 
we replace 34 by a smaller set, of c n — z k>2n — Vk,2n elements, which together with the 
order In — 1 elements of Z k , forms a basis of the c n elementary differentials. In terms of 
the generating functions, we add the final step to the iteration of Eq. (|l 

Vk+i 



^2 min(l/fc+i,2n> c n ~ z k +i,2n)t 2n (20) 



n>l 



Let the resulting limiting formal series be x E (t), VE(t), and z E (t). The generating function 
for XE(t) is then computed to be 

x E (t) =t 2 + t A + f + 2t 8 + 3t 10 + 6t 12 + lit 14 + 23t 16 + 47t 18 + 102t 20 + 221t 22 + 

484t 24 + 1069t 26 + 2386t 28 + 5364t 30 + 12143 32 + 27645t 34 + 63259t 36 + . . . 

(21) 

which should be compared with Eq. (p"4]). At order 14, 16, and 18 the dimensions are 
limited by the number of elementary differentials, but for n > 9, [t 2n ]xE{t) < [t 2n ]x(t) < 
c n . Because the new map on generating functions, Eq. (|2"0|), is not analytic, it is harder to 
determine the location of its smallest singularity. We found the smallest root of successive 
polynomial truncations of 1 — t — and extrapolated these results to obtain 

1/^ = 1.8250339..., 1/^ = 1.574.... (22) 

These are upper bounds for the entropy of the class of Euclidean mechanical systems and 
their modified potentials, respectively. 

(Murua |2(J has also considered this case, in the context of order conditions for Hamil- 
tonians of the form \ ^p 2 + V(q). He finds a unique independent tree of a certain type 
for each order condition, and enumerates these up to order 6. It would be interesting to 
compare the two approaches at higher order.) 

The situation is quite different for non-Euclidean, i.e. general, mechanical systems. 
Repeating the above calculation for a general kinetic energy T(p) = ^p T M(q)p, we get 
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the following modified potentials. The associated trees will be explained below. 
Z x = V = • 

Z 2 = [BBA] = M(V', V) = S\ 

Z 3 = [BBA, BA] = 2M(V, V"{M{V'))) + M'(V, V, M(V')) = 2 JV\ + 
Z 4 = [BBA, [BBA, A]] 

= AV'(M(V"(M(V"(M(V'))))) + 3V'(M(V"(M(M'(V, V')))))+ 
M{M'{V',V'),M'{V',V')) 

= jy^ + ><>< (23) 

Z 5 = [[BBA,BA},BA] 

= AV'(M(V"(M(V"(M(V')))))) + 2V'"(M(V), M(V), M(V'))+ 
6M'(V, M(V), V"{M{V'))) + M'(V, V, M(V"(M(V'))))+ 
M'iV', M(V), M'(V, V')) + M"(M(V'), M(V), V , V) 

In this case each modified potential of order 2n — 1 is a scalar elementary differential of 
V and M. These correspond to bicoloured free trees with 2n — 1 nodes, of which n nodes 
are labelled V (shown as solid circles above) and n — 1 nodes are labelled M (shown as 
open circles above); the latter must have at least 2 branches since a derivative of M has 
at least 2 indices. Of the 1, 1, 3, 11, 47, and 235 free trees of order 1, 3, 5, 7, 9, and 11 
respectively, exactly 1, 1, 2, 8, 34, and 175 of them can be coloured (labelled) in this way. 
The calculation above shows that of these colourable trees, precisely one colouring of 1, 
1, 2, and 7 of these colourable trees occur in the modified potentials of orders < 7. (The 
other colourings of these trees do not occur, because of the way in which the trees at each 
order are built from the trees of lower order. The colourable 7-node tree also does 

not occur.) It is clear that there is enormously much more freedom in this case than in 
the ('Euclidean', T(p) = |E 

pf) case considered previously. Therefore, we believe that 
all the modified potentials are independent in this case. This supports Conjecture |Tj 

5.3 Other polynomially graded Lie algebra 

We close with a list of some other Lie algebras of class ^3. In each case one can consider 
the case of two generators A and B of degrees 2 and and the induced homomorphism 
from Ly(A, B). 

1. The case of classical mechanics. The objects are real functions on a cotangent 
bundle, homogeneous polynomial in p. This can be specialized to the following 

cases. 
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(a) Q any Riemannian manifold, any potential energy, degree(X) is the total degree 
of X in p. Entropy is < 1.8254..., Eq. (|18D , with Conjecture p] implying 
equality. 



(b) Q = R n with the Euclidean metric. Entropy is < 1.8250. . ., Eq. (|22|). It is 
remarkable that these two Lie algebras, not previously distinguished from each 
other in the literature, differ starting at order 13, and have slightly different 
entropy. 

(c) Q = ~R n , functions polynomial in p and q. We can then introduce a bigrading by 
degree in p and by degree in q. To get a new Lie algebra, one of the generators 
has to be degree in each grading, which forces Q Euclidean, A = iJ^Pi, 
B = V(q) polynomial. For example, we have computed the dimensions of 
the Lie algebra generated by cubic potentials for small n in Table 2 — they are 



remarkably small. See |l 1] for an analysis of this case in terms of special types 
of trees. 

Homogeneous polynomial vector fields on R m graded by total degree mxi,...,xj t for 
some 1 < k < m. In the case k = m, the vector fields in X associated with L(A, B) 
(degree(v4) = 2, degree(-B) = 0) are associated with free trees in which each node 
has degree at most 2 (since only the first two derivatives of A are nonzero). Their 
numbers are 1, 1, 1, 2, 3, 6, 11, 23 (so far the same as for the free trees), then 



46, 98, 207, 451, .... (Sloane's A001190 p2[), which gives an upper bound for the 



number of independent elements of Z of each odd order. These grow more slowly 
than the free trees, and even more slowly than 2, with entropy 1.5758, compared to 
1.5821 (Eq. (li~9])) for Z. Perhaps in this case the trees X generate the Lie algebra 
as 1®L{A, [1, A})1 

3. As the previous item, but multigrading by total degree in different subsets of the 
variables. 

4. Homogeneous polynomial vector fields with the variables partitioned (x, y) with 
x E y E ]R m , and the vector fields of the form p-^ + q-^ with either degree y (g) < 
degree y (p) + 1, or p = and degree J/ (g) = 0. Simple mechanical systems form 
examples of this class. So do high-order ODEs of the form y^ = f(y, . . . ,y( n ~ 2 ^) 
when re-written as first-order systems 

ii = Xi+i, i = 0, . . . , n — 2, 
i n -i — f(%i, ■ ■ ■ > 2^-2); 

with Xi = y^\ k = n — 1, and m = 1. 



21 



5. Consider the Schrddinger equation 

itjj = vV> + V(x)ip, 

where V 2 is the Euclidean Laplacian. The two operators V 2 and V(x) generate a 
Lie algebra of class where the grading is by degree of the differential operators. 
For example, 

[v 2 ,u]^ = v- (vi>) + vv- i> 

is of degree 1, 

[U,U,V 2 h/> = (V • (V 2 ))^ 

is of degree 0, and 

[V,V,V,V 2 }^ = 0. 
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